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W i ABSTRACT 

<N . 

We have used the Arecibo telescope to measure the HI absorption spectra of 
eight pulsars. We show how kinematic distance measurements depend upon the 
values of the galactic constants R Q and O , and we select our preferred current val- 



ues from the literature. We then derive kinematic distances for the low-latitude 
pulsars in our sample and electron densities along their lines of sight. We com- 
bine these measurements with all others in the inner galactic plane visible from 
Arecibo to study the electron density in this region. The electron density in the 
interarm range 48° < I < 70° is [0.017 (-0.007, +0.012) (68% c l)} cm~ 3 . This 
i s 0.75 (-0.22, +0.49) (68% cl) of the value calculated by the ICordes fc Lazio 



( 120021 ) galactic electron density model. The model agrees more closely with elec- 



OO 

m 

^ \ tron density measurements toward Arecibo pulsars lying closer to the galactic 

center, at 30° < / < 48°. Our analysis leads to the best current estimate of the 
distance of the relativistic binary pulsar B1913+16: d = (9.0 ± 3) kpc. 

We use the high-latitude pulsars to search for small-scale structure in the in- 
terstellar hydrogen observed in absorption over multiple epochs. PSR B0301+19 



exhibited significant changes in its absorption spectrum over 22 yr, indicating HI 
structure on a ~ 500 AU scale. 

Subject headings: Pulsars: distances — Galaxy: fundamental parameters — ISM: 
structure — ISM: abundances — ISM: clouds — Radio lines: ISM 
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Introduction 



Neutral hydrogen (HI) absorption measurements of pulsar signals at A = 21 cm are 
important probes of va rious properties of the interstellar medium s uch as the small-scale 
structure of of cold HI (jDickey et al.lll98ll ; IStanimirovic et al.ll2003bl ) and c alibrators of the 



pulsar distance scale and electron density models at large galactic distances (jWeisberg et al. 



19791 . 119871 . Il995l ; iFrail fc Weisbergill99Cll ). These measurements are compleme ntary to inter 



ferometrically determined pa rallaxes, which can be utilized on nearer sources (IBrisken et al. 



20021 ; IChatterjee et al.l 12004 ) . We report new Arecibo HI absorption observations of eight 
pulsars, and we use these measurements to study the electron density of the interstellar 
medium and the small-scale structure of neutral hydrogen clouds. 

Sensitivity limitations indicate that few if any additional pulsar HI absorption measure- 
ments of electron density in the galactic plane at the intermediate first-quadrant galactic 
longitudes accessible to the Arecibo telescope will be made in the next decade or so. Hence 
it is an appropriate time to combine the new results with all previous pulsar HI measurements 
at these longitudes to globally assess the density in the galactic plane in this region. We will 
use this information to estimate the distance to the relativistic binary pulsar B1913+16. In 
addition, HI absorption measurements on some of our high-latitude pulsars were originally 
measured over twenty years ago, during which interval the pulsars have moved many AU 
through the interstellar medium. Hence comparison of the old and new absorption spectra 
yields information on the small scale structure of the absorbing neutral hydrogen. 

The plan of the paper is as follows: The pulsar HI absorption observational technique is 
sketched in the next section (§2). We present HI absorption spectra and kinematic distance 
results for low-latitude pulsars in §3. In §4, we use all measured pulsar distances in the 
inner galactic plane visible from Arecibo Observatory to analyze the electron density in that 
region. In §4.1, we review the latest work on the galactocentric distance of the Sun and the 
orbital velocity of the Local Standard of Rest in order to select an optimum model to use in 
refining pulsar kinematic distances and galactic electron densities. Then in §4.2, we apply 
the results of §4 and §4.1 to determine the distance of the relativistic binary pulsar PSR 
B1913+16. We provide absorption spectra of high-latitude pulsars in §5, along with analyses 
of time- variability of absorption in those cases where earlier epoch data are available. Finally, 
we discuss our conclusions in 86. 
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Observations 



All observations were made with the 305-m Arecibo telescope from 1998 to 2000. The ra- 
dio frequency signals near 142 MHz were mixed to bas eband, sampled, and recorded with the 
Caltech Baseband Recorder (ICBR. Jenet et al.ll200ll ). a 10 MHz bandwidth, fast-sampling 
receiver backend. Every 100 ns, the CBR sampled complex voltages with four-level digitiza- 
tion from the two orthogonally polarized feed channels, and recorded the samples to tape for 
later processing. (For additional details of the observing techniques and e quipment, see Sta- 
nimir ovic et al [2003a,b]). The data were then corrected for quantization (IJenet fc Anderson 
19981 ). Fourier transformed and folded modulo the apparent pulsar period with the Super- 
computers of Caltech's Center for Advanced Computing Research, resulting in data cubes 
consisting of intensity as a function of pulsar rotational phase (128 phase bins) and radio 
frequency (2048 frequency channels, each having 1 km/s width). Subsequent processing at 
Carleton College collapsed the data cubes into two spectra: the pulsar- on spectrum is a sum 
of the spectra gathered during the pulsar pulse weighted by I psr ((f)) 2 , while the sum of those 
gathered in the interval between pulses is called the pulsar- off spectrum. Here I psr {4>) is the 
broadband pulsar intensity in a given pulse phase bin <fi. 

Two final spectra are formed and displayed for each pulsar. The pulsar absorption 
spectrum, which represents the spectrum of the pulsar alone less any absorption caused by 
intervening neutral hydrogen, is the normalized difference of the pulsar- on and -off spectra. 
In order to maximize sensitivity, multiple integrations are summed with a weight depend- 
ing on the square of the pulsar signal strength I psr (t) during each integration. In some 
cases, additional sensitivity was achieved by Hanning smoothing the final absorption spec- 
trum, yielding a resolution of 2 km/s. (Any absorption spectrum that has been Hanning 
smoothed is labelled as such when displayed below.) The HI emission spectrum is the time- 
integrated, pulsar- off spectrum, calibrated in bright ness temperature units by matching its 



peak with the Leiden /Argentine /Bonn HI S urvey (IHartmann fc Burton! 119971 ; lArnal et al. 



2000l ; iBajaja et al.ll2005l ; iKalberla et al.ll2005l ). All spectra were frequency switched and low- 
order polynomials were fitted to and removed from them in some cases in order to flatten 
intrinsic or scintillation-induced bandpass ripples. 

The basic observing parameters are tabulated for low- and high-latitude pulsars in the 
early columns of Tables 1 and 3, respectively. The quantities T sysfi{i _\ ine and ov^ff-iine, the 
system temperature and measured 1-a noise in optical depth away from the HI line, are 
both given. The former value is the sum of an estimated 40 K receiver contri bution, plus a 



sky b ackground determined by extrapolating the 408 MHz sky temperature (IHaslam et al. 



19821 ) at the pulsar position to 1420 MHz with a spectral index of -2.6. The HI emission line 
itself contributes significantly to the system temperature (and hence to the noise) at those 
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velocities where it is present. To determine the expected optical depth noise at any velocity 
v, one may use cr T (v) = cr ri0ff _ linc x [T HI (v) + T sySi0ff „ line ]/T sySi0ff _ line , where T HI (v) is the 
measured brightness temperature of HI at that velocity. 



3. Kinematic Distance Analyses and Results on Low-Latitude Pulsars 



A low-latitude pulsar HI absorption spectrum can be used to set limits on the pulsar's 
distance, using the "kinematic" technique. A pulsar is farther than an HI cloud that absorbs 
its signal, and closer than one that does not. The latter "no absorption" criterion is impossi- 
ble to ascertain in real spectra because of the inevitable presence of noise which could mask 
weak absorption. However it was shown by IWeisberg et al.l (119791 ) that emission features 
with Tb > 35 K almost always exhibit significant absorption of radiation from background 
objects. Therefore subsequent investigators have assigned an upper distance limit only at the 
velocity where o ne finds both no significant absorption, and an emission feature of TJ, > 35 
K. We use a flat iFich et al.l (119891 ) galactic rotation model an d the IAU galactic constants of 
R = 8.5 kpc and G D = 220 km/s (IKerr &: Lynden-Belllll986l ) to convert the radial velocities 
to distance. The resulting model radial velocities are shown in a panel under each absorp- 
tion/emission spectrum pair. We also add and subtract velocities of 7 km/s to our nominal 
velocities, in order to derive estimates of the uncertainties in distance limits due to stream- 
ing and random gas motions in the Galaxy. These procedures a re identical to tho s e use d 
in the critical evaluation of all such measurements then extant by lFrail fc Weisbergi (ll990l ). 
and by all subsequent pulsar HI absorption experimenters. Hence our results are directly 
intercomparable with these earlier ones. Our derived distance limits are discussed below for 
each pulsar, and are summarized in Table 1. See §4.1 for a discussion of modifications to 
our derived distances for different values of the Galactic constants. 



3.1. PSR J1909+0254=B1907+02 (/=37?6; b=-2°7), Fig. Q] 



The farthest observed absorption feature is at v — 60 km/s, well before the tangent 
point. Hence the pulsar lies beyond a lower distance limit of 3.8 ± 0.5 kpc. Unfortunately 
the HI emission has T& < 10 K at all velocities where one might test for an upper distance 
limit via lack of absorption, because the b = 2?7 line of sight rapidly leaves the hydrogen 
layer. Since signific ant absorption could not in any case be guaranteed at the velocity of 
such weak emission (IWeisberg et al.lll979l ). no upper distance limit can be set. 
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3.2. PSR J1922+2110=B1920+21 (/=55?3; 6=2?9), Fig. H 



Weisberg et al.l (119871 ) observed this pulsar. Their HI spectrum was contaminated by 



interstellar scintillation and was rather noisy, pr ompting us to r eobser ve it with greater 
sensitivity. The highest velocity absorption that IWeisberg et al.l (119871 ) could confidently 
detect was at v = 26 km/s, leading to their limit of (1.9 ± 0.7 kpc ) < d. Our new 
observations clearly exhibit absorption at v = 41 km/s, near the tangent p oint. Hence we 



revise the lower distance limit significantly upwards: (4. 8 ±1.8 kpc ) < d. The lWeisberg et al. 



( 119871 ) spectrum also showed an absorption feature at this velocity but it was not sufficiently 
above the noise to serve as a reliable kinematic distance indicator. The dip in the current 
absorption spectrum at v — —48 km/s is probably noise and will not be used to set a distance 
limit. (It appears much more prominently in one of the two c ircula r polarization channels 
than in the other, and is not visible in the IWeisberg et al.l (119871 ) absorption spectrum. 
Furthermore, the nearby source G55.6+2.3=B1923+210 does not exhibit reliable far side 
absorption until v < —60 km/s [Dickey et al 1983; Colgan et al 1988].) Conversely, the lack 
of absorption in the TJ, = 41 K emission feature at v — —65 km/s sets the upper distance 
limit of d < (16.2 ± 1.0) pc. 



3.3. PSR J1926+1648=B1924+16 (/=51?9; 6=0?1), Fig. E3 



While the Fich et al. model predicts a tangent point velocity of 48 km/ s in this direction, 
we detect significant emission and absorption well beyond this velocity, with the las t major 
feature c entered near v = 61 km/ s. T hese features are due to the Sagittarius arm (IBurton 



1970J). iGarwood fc Dickeyl (119891 ) and IColgan et al.l (119881 ) also observe HI absorption out 
to these velocities in the nearby sources G50. 625-0. 031 and G51. 4-0.0, respectively. We 
choose the tangent point as our lower distance limit, so (5.2 ± 1.8) kpc < d. The lack of 
absorption in the TJ, ~ 42 K emission feature at v ~ —47 km/s yields the upper distance 
limit: d < 14.9 ± 0.8 kpc 



4. Analysis of the Electron Density in the Inner Galactic Plane Visible from 

Arecibo 

The currently reported HI absorption distance measurements are probably among the 
last to be determined at Arecibo in the foreseeable future, as all known pulsars having suf- 
ficient flux density to achieve adequate S/N in a reasonable time have now been measured 
as well as is possible with this procedure. Unfortunately galactic HI emission is now the 
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dominant noise source, so that future receiver improvements will not significantly decrease 
the overall noise. Therefore, new low-latitude pulsars discovered in the future in this lon- 
gitude range will probably be too faint for the HI absorption technique, even at Arecibo. 
Consequently now is a good time to review all such measurements and their implications for 
the electron density in the inner galactic plane visible from Arecibo. 

The dispersion measure DM of a pulsar, derived from multiple-frequency timing mea- 
surements and reported in Manchester et al. ( 2005 1 "*! directly yields the path- integrated 



electron density along the line of sight: DM = J n e ds. Hence our distance limits, d, can be 
combined with the published dispersion measures to yield mean electron densities (n e ) along 
the line of sight: 

Jn e ds DM 

Table 2 lists the electron density limits derived in this fashion from our measurements 
and from all earlier ones in similar directions; i.e., toward the inner galactic plane accessible 
to the Arecibo telescope. Note that the kinematic distances listed in Table 2 were all derived 
wit h the standard pr o cedur es discussed above in §3, namely the unifor m criteria established 



by iFrail fc Weisberd (119901 ) and the IAU standard galactic constants (IKerr fc Lynden-Bell 



19861 ). Procedures to modify these values for other choices of galactic constants are given 



below in §4.1. 

In Fig. Hk, we display measured limits on the mean electron density as a function of 
galactic longitude I in the inner galactic plane visible from Arecibo. Only those pulsars from 
Table 2 with \b\ < 3° and a lower distance limit d\ > 1 kpc are included in Fig. 0Ji in order 
to concentrate on kpc-scale averages in the galactic plane. Almost all of these measurements 
were derived from Arecibo pulsar HI absorption spectra. While there is significant scatter, 
it is apparent that the measured densities tend to decline as the line of sight rotates away 
from the inner Galaxy, with a notable drop in the region between the Sagittarius arm (the 



first a rm interior to our location) and the local spiral arm (see Fig. [5]). lAbles fc Manchester 



(119761 ) were the first to have adequate (n e ) measurements to infer that interarm densities 
are lower than those in spiral arms. As illustrated in Fig. Hk, the numerous measurements 
made since that epoch serve to confirm and refine their suggestion, at least in this particular 
interarm region. Indeed, the drop in (n e ) as I exceeds 48° corresponds with the line of 
sight finally reaching a longitude where it no longer intercepts the Sagittarius arm (see Fig. 
[5]). Meanwhile, with only one exception, every pulsar measurement in the interarm range 



1 The catalog is maintained and updated at http://www.atnf.csiro.au/research/pulsar/psrcat 
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48° < / < 70° is consistent with a relatively low (n e ) ~ 0.02 cm -3 (see below for further 
analysis). The electron density at lower latitudes sampled in Fig. is clearly significantly 
higher, with an average in the vicinity of 0.05 cm -3 (plus superposed variations) in the 
30 — 40° longitude range. 

It is also useful to compare our updated set of electron density measurements with the 
most widely used galactic density model (NE2001, Cordes & Lazio 2002). Fig. Hb exhibits 
the measured limits on (n e ), normalized by the mean electron densities predicted by the 
NE2001 model. The measured-to-model electron density ratio hovers near 1 (with some 
inevitable scatter) for I < 48°, indicating that the model fits the data adequately at those 
longitudes. However, at the higher longitudes 48° < / < 70° discussed above, there is a 
tendency for the typical measured-to- model density ratio to lie below 1, suggesting that the 
NE2001 model densities are somewhat high in this interarm region. 

In order to further explore the electron density in this interarm region, we performed 
Monte Carlo simulations to assess the best value and uncertainty on (n e ) and on the 
measured-to-model (n e ) ratio for the pulsars in this region. Of the ten pulsars with measured 
distance limits and limits on (n e ) in this range (see Table 2), we discarded the data from PSR 
J1935+1616 since it lacks a measured upper distance and lower density limit. For each of the 
remaining nine pulsars, all possible distances between the measured upper and lower limits 
were rendered equally probable in our simulations by choosing distances randomly from a uni- 
form distribution between the measured distance limits. Note that the assumption of uniform 
probability is the simplest reasonable hypothesis for measurements such as these which have 
only a lower and upper limit. We find that (n e ) = [0.017 (-0.007, +0.012) (68% c.l.)] cm" 3 , 
and that the ratio of measured-to-NE2001 electron density is 0.75 (-0.22, +0.49) (68% c.l.) 
for 48° < I < 70°. Future models should take this discrepancy into account. 

For simplicity, both the measured and model densities discussed above are based on the 
standard IAU galactic constants. As described below in §4.1, the ratio plotted in Fig. Hb 
remains virtually identical if both constituents are rescaled to incorporate other choices of 
galactic constants. 

4.1. The Effects of Improved Galactic Constants R a and O 

Much progress has been made on galactic structure and kinematics in the ~ 20 years 
since the IAU constants were defined. Therefore it is useful to discuss improved measure- 
ments of R and and 9 Q , and to investigate the influence of these better values on the 
kinematic distances and electron densities discussed in §3 and §4. 
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From remarkable proper motion and radial velo city measurement s of th e orbit of star 



isenhauer et ali (2005h fi nd that the 



S2 about Sgr A* through much of its 15-yr period 

distance to the galactic center R Q = 7.62 ± 0.32 kpc. iReid &: Brunthalerl (l2004j ) have made 
equally stunning interfero metric measurements of the proper motion of Sgr A*, which yield 
an angular velocity of the Local Standard of Rest (LSR) about the galactic center, Q Q = 
o /R o = (29.45 ± 0.15) km s _1 kpc -1 . Both of these measurements are notable in that 
they rest on far fewer assumptions than earlier determinations, and hence should be freer of 
systematic errors. Consequently we adopt them in what follows. An even newer estimate of 



R n fr om infrared measurements of red clump giants in the galactic bulge (INishiyama et al. 



20061 ) yields a value of 7.52 ± 0.10 (statistical) ±0.35 (systematic) kpc, which is consistent 
with our above adopted value. 

Consider material traversing a circular orbit in the plane at galactocentric radius R. 
Given a flat rotation curve with circular velocity O , its radial velocity with respect to the 
LSR is v ra d = <d {Ro/R — l)sin/. With R = [R 2 + d 2 — 2R dcosl) 1 ^ 2 , the distance as a 
function of v ra d is given by 



d = Ro( cos/ ± Jl + - - - sin 2 / . (2) 
\ V B sm/ / 

Hence our new adopted galactic parameters, giving a small change in o but a relatively 
large change in R , result in a recalibration of the distance in a particular direction, accom- 
plished by rescaling it with an approximately constant factor of R ,new/Ro,oid = 7.62/8.50 = 
0.896. Kinematic distance limits and galactic distance models that use the old IAU galactic 
constants should be adjusted by this factor to reflect the improved constants. 

Fig. [6] displays distance - radial velocity curves toward / = 50° derived with the old and 
new flat galactic rotation curves. Fig. [6^, illustrates that the old and new model distances 
differ primarily by a constant multiplicative factor, as indicated in the previous analysis. 
It follows from the above, that when distance is normalized by R Q , the old and new radial 
velocity - distance curves are almost identical (Fig. [6b). Such normalized distances are 
useful in the analysis of the acceleration of pulsars in the galactic gravitational field, as will 
be discussed in the next section (§4.2). 

As described above, the new galactic constants would result primarily in the measured 
distances being rescaled by the ratio of new to old R Q . The electron densities, including model 
densities that have been calibrated via pulsar kinematic distance measurements, would then 
be rescaled by the inverse of this ratio (see Eq. 1). Hence Fig. would be thus rescaled; 
while Fig. Hb would remain unchanged since it displays the ratio of two densities, each of 
which should be rescaled by the same factor. 
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4.2. The Distance to the Relativistic Binary Pulsar B1913+16 



The relativistic binary pulsar B1913+16 lies in the heart of the region of the Galaxy that 
we are studying, at (I, b) ~ (50°, 2°). The orbital decay of this s ystem due to gravitationa l 
wave emission provides a strong test of relativistic gravitation (ITaylor fc Weisberd Il989l ). 
Currently, the precision of the relativistic tes t is limited by uncerta inties in the galactic 
acceleration of the pulsar and the solar system (IDamour fc Taylorlll99ll ). which are currently 
dominated by uncertainties in t he pulsar distance d [or more s pecifically, uncertainties in S, 
the distance normalized by R Q (jWeisberg fc Taylor! 120031 . 120051 )]. Hence it is very important 
to assess the current results in an effort to improve our knowledge of B1913+16's distance. 

Fig. shows that the measured electron density appears to have a local minimum near 
the pulsar's longitude. The relatively low value of (n e ) at / ~ 50° is not surprising, since we 
noted in §4 that the line of sight in this direction traverses a long path between spiral arms (see 
Fig. The tightest limits on (n e ) in the whole inner plane visible from Arecibo are for PSR 
B1915+13 at I = 48°, which happens also to be the pulsar closest to PSR B1913+16. Hence 
we adopt the measured limits from B1915+13, including an allowance for up to 7 km s" 
circular velocities (not shown in Table 2 or Fig. H]): (n e ) ~ (0.0188 ± 0.006) 



..5 kpc 
Ro 



non- 

3 



cm 



Note that this value is congruent with and has modestly tighter limits than the electon 
density derived in this interarm region in §4. 

Armed with our new estimate of (n e ), we use Equation 1 to find the distance to PSR 
B1913+16, which has DM ~ 169 pc cm" 3 : 



<W 3+ i6 = (10.0 ± 3.2) 



R 



8.5 kpc 



(3) 



kpc. Our result can also be compared with the ICordes fc Laziol (120021 ) NE2001 model dis 
tance of 5.90 ( ^,°„ „ ^ kpc, which is lower for the reasons discussed above. Damour fc Taylor 

ftp 



.5 kpc 



( 1991 ) used techniques similar to ours to estimate that G? B i9i3+i6 = (9-2 ± 1.4) ^ 



.5 kpc 

While the error bars of their and our esti mates overlap sign i ficantl y, we are not able to 
justify tightening them to the degree that IDamour &: Taylor! (Il99ll ) did. Finally, adopt- 
ing what we judged in §4.1 to be the current best measurement of R (R Q = 7.62 ± 0.32 
kpc; Eisenhauer et al. 2005) leads to our best estimate of the distance to PSR B1913+16: 
4si9i 3 +i6 = (9.0 ± 3) kpc. 

As noted above, the uncertainty in the general relativistic orbital decay rate of this pulsar 
is currently dominated by the uncertainty in the quantity #Bi9i3+i6- For these purposes, it 
is useful to give 

5 B i9i3+i6 = (d B1913+16 /R ) = (1.18 ± 0.38). (4) 
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High Galactic Latitude Pulsars and Searches for Temporal Variability of 

Absorption 

Relative ly nearby pul sars are excellent targets for studies of the tiny-scale atomic struc- 



ture (TSAS, iHeilesI ( 119971 )). Comparison of their HI absorption spectra at multiple epochs 
as they move across the sky permits us to study structure down to AU scales. We discuss 
here our measurements of five such pulsars. All but two have been previously observed, so 
we will search for changes between the previous and current epochs. The pulsars' observing 
parameters, distances, and transverse velocities are listed in Table 3. 

The pulsars discussed in this section are all located at high galactic latitudes (in all 
cases, |6| > 25°). Since their lines of sight leave the galactic hydrogen layer fairly quickly, 
they are not amenable to the kinematic distance technique discussed in §3. Conversely, we 
can not search for temporal absorption variations in the low-latitude pulsars of §3 because 
two have never before been observed, and the third's absorption spectrum was too noisy 
at the earlier epoch to yield meaningful differences. (Indeed we reobserved it in order to 
more confidently measure its kinematic distance because the noise in the earlier spectrum 
precluded a secure result.) 



5.1. J0304+1932=B0301+19 (/=161?1; 6=-33?3), Figs. H and El 



The HI absorption spectrum of this pulsar was first observed in 1976-77 by lDickey et al. 



(11981I ). That spectrum and our new one are plotted together in Fig. [71 Both epochs' spectra 
have a similar ~ 1 km/s resolution. Both the integrated optical depth (the equivalent width) 
and the depth of most individual absorption components have changed significantly between 
the two epochs. Although the early observations were gathered with a spectrometer that 
digitized the signal t o only one bit, i ts response to rapidly varying pulsar signals was very 
well-characterized by Weisberj (19781). Careful q uantization corrections have also been made 
for the current epoch (IJenet fc Anderson! Il998l ). Hence we believe that both spectra are 
accurate and that the change is real. 

The channel spacing of the old and new spectra are 1.056 and 1.031 km/s, respectively. 
In order to study the absorption changes more carefully, we slightly coarsened the spacing of 
the new data via interpolation to match the channel spacing and center frequency of the old 
data. We also ensured that the old and new channel frequencies matched by verifying that 
the (essentially noise-free) emission spectra were consistent to less than a channel width. 
Fig. [S] displays the old and (resampled) new absorption spectra (top) and their difference 
(bottom). The difference spectrum also displays a ±2cr noise envelope in order to assess 
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the significance of variations. [The noise envelope grows at the cent ral velocities because 



the HI emission line itself contributes significantly to the system noise ( IJohnston et al.ll2003 
Stanimirovic et al.ll2003bl ).] The difference spectrum clearly shows a general trend whereby 



the absorption line depth is greater at the later epoch in most of the central channels, with 
the largest single- channel difference being significant at the 2.6a level. 

The ~ 22 year time baseline is unique in being one of the longest extant showing 
absorption variations, leading to a length scale of ~ 500 AU which is second only to the PSR 
1557-50 scale of ~ 1000 AU (IJohnston et al.l 120031 ) . Further implications of these results 



for small-scale structure in the interstellar medium will be analyzed in a separate paper 
(Stanimirovic et al, in preparation). 



5.2. J1239+2453=B1237+25 (/=252?5; 6=86?5), Fig. H 



Dickey et al.l (1198ll ) also measured the HI absorption spectrum of this pulsar in 1976-77. 



The original and current pulsar spectra show no absorption. These results are not surprising, 
given the galactic polar line of sight and the very weak HI emission in this direction. The 
pulsar moved ~ 2200 AU between the two observations. 



5.3. J1537+1155=B1534+12 (/=19?8; 6=48?3), Fig. [TO] 

This pulsar is a mem ber of a double ne utron star binary system. Its pulse timing par- 
allax gives d > 0.67 kpc (jStairs et al.ll2002l ). Under the assumption that general relativity 
provides the correct description of gravitational wave emission, the excess observed orbital 
period change not attributable to gravitational wave damping yields d = (1.02 ± 0.05) kpc 
( jStairs et al.ll2002l ). There is no statistically significant absorption evident in our pulsar spec- 
trum along this high galactic latitude line of sight. No previous HI absorption measurements 
have been made on this source. 



5.4. J1543+0929=B1541+09 (/=17?8; 6=45?8), Fig. UJ 

No statistically significant absorption is present in this high galactic latitude pulsar's 
spectrum. This is the first absorption spectrum from this pulsar. 
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5.5. J2305+3100=B2303+30 (/=97?7; 6=-26?7), Fig. [H 



The broad, shallow dip in the displayed spectrum disappears in one of our two orthogonal 
polarizatio ns, leading us to co nclude that it is not real. This pulsar was also observed in 
1976-77 by lDickey et al.l (ll98ll ). In that case as well, no significant absorption was seen. The 
pulsar traversed ~ 1600 AU in the intervening time. 



6. Conclusions 

We have determined the HI absorption spectra of eight pulsars. The three low-latitude 
pulsars yield kinematic distances and electron densities in the inner galactic plane visible from 
Arecibo. These observations mark the completion of a two-decade effort to accurately mea- 
sure the HI absorption spectrum of all such pulsars that are strong enough to be accessible 
to this technique. Therefore, we have combined our new measurements with all others in this 
direction to study the electron density in this region. We find that the mean electron density 
in the plane in the interarm range 48° < / < 7 0° is [0.017 (-0.007, +0 .012) (68% c.l.)} cm" 3 , 



which is 0.75 (-0.22, +0.49) (68% c.l.) of the ICordes fc Laziol (120021 ) model value currently 



used by most re s earch ers. At the lower longitudes accessible to Arecibo (30° < I < 48°), the 



Cordes fc Laziol (120021 ) model appears to conform generally to the measurements, aside from 



expected local variations. 

As part of the process, we show how to modify kinematic distances and electron densities 
as a function of the galactic constants R Q and O . We review recent efforts to determine 
the values of these constants and select the best current answers. Applying all of these 
results to the relativistic binary pulsar B1913+16, we find a dispersion measure distance of 
d = (10.0 ± 3.2) ( s.ifi p P ) k pc, or d = (9.0 ± 3) kpc if we adopt our current best choice for 



R (jEisenhauer et al.l 120051 ) 



The five high-latitude pulsars are most useful in multiepoch studies of small-scale struc- 
ture in the interstellar medium. Two were observed for the first time by us in this experiment. 
Of the other three, two showed no measurable absorption at either of two epochs separated 
by twenty-two years; while one, PSR B0301+19, exhibited a significant change in absorption 
profile over that timespan, indicating HI structure on a ~ 500 AU scale. 



JMW, AH, and AF were supported by grant AST 0406832 from the National Science 
Foundation. Arecibo Observatory is operated by Cornell University under cooperative agree- 
ment with the NSF. 
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Table 1: Observing Parameters and Measured Kinematic Distances for Low Latitude Pulsars 



PSR J 


PSR B 


^obs 


^~sj/s, of? — line 


°Y,off — lino 


I 


b 


DM 














(hr) 


(K) 




(dog) 


(deg) 


(pc cm -3 ) 


(km/s) 


(km/s) 


(kpc) 


(kpc) 


J1909+0254 


B1907+02 


3.0 


45 


0.11 


37.6 


-2.7 


171.7 


60 




3.8 ±0.5 




J1922+2110 


B1920+21 


2.5 


43 


0.09 


55.3 


2.9 


217.1 


Tpc 


-65 


4.8 ± 1.8 


16.2 ± 1.0 


J1926+1648 


B1924+16 


7.4 


45 


0.06 


51.9 


0.1 


176.9 


T pc 


-47 


5.2 ± 1.8 


14.9 ±0.8 



a vi and v u are lower and upper velocity limits, respectively. 

b di and d u are lower and upper dista nce limits, respectively. A flat lFich et alj ()1989f ) rotation curve with the 



standard IAU values of R and Q a |Kerr fc Lvnden-Bell 19861 ) was used. 
c The tangent point velocity is adopted as the lower limit. 



Table 2: Measured Electron Densities toward Pulsars with 30° < / < 70°; \b\ < 3°; and 1 
kpc < d. 



PSR J 


PSR B 




I 

(deg) 


b 

(deg) 


DM 

(pc cm -3 ) 


dl a 
(kpc) 


(kpc) 


; — 

(n e , u ) b 

(cm -3 ) 


; rn — 

(cm -3 ) 


Distance 
Ref. c 


Assoc. 


J1848-0123 


B1845-01 


31.34 


0.04 


159.5 


4.2 


4.8 


0.038 


0.033 


FW90 




J1852+0031 


B18494 


-00 


33.52 


0.02 


787.0 


7.1 


16.6 


0.111 


0.047 


FW90 




J1856+0113 


B18534 


-01 


34.56 


-0.50 


96.7 


2.7 


3.9 


0.036 


0.025 


C75 


SNR:W44 


J1857+0212 


B18554 


-02 


35.62 


-0.39 


506.8 


6.9 




0.073 




FW90 




J1901+0331 


B18594 


-03 


37.21 


-0.64 


402.1 


6.8 


15.1 


0.059 


0.027 


FW90 




J1901+0716 


B18594 


-07 


40.57 


1.06 


252.8 


2.8 


4.7 


0.090 


0.054 


FW90 




J1902+0556 


B19004 


-05 


39.50 


0.21 


177.5 


3.1 


4.3 


0.057 


0.041 


FW90 




J1902+0615 


B19004 


-06 


39.81 


0.34 


502.9 


6.5 


15.8 


0.077 


0.032 


FW90 




J1903+0135 


B19004 


-01 


35.73 


-1.96 


245.2 


2.8 


4.0 


0.088 


0.061 


FW90 




J1905+0154A 






36.20 


-2.20 


194.0 


6.4 


8.2 


0.030 


0.024 


R97 


GC:NGC6749 


J1906+0641 


B19044 


-06 


40.60 


-0.30 


472.8 


6.5 


14.0 


0.073 


0.034 


FW90 




J1909+0254 


B1907+02 


37.60 


-2.71 


171.7 


3.8 




0.045 




this work 




J1909+1102 


B19074 


-10 


44.83 


0.99 


150.0 


1.3 


6.0 


0.035 


0.025 


FW90 




J1915+1009 


B1913-) 


-10 


44.71 


-0.65 


241.7 


6.0 


14.5 


0.040 


0.017 


FW90 




J1916+1312 


B1914-( 


-13 


47.58 


0.45 


237.0 


1.0 


5.7 


0.059 


0.042 


FW90 




J1917+1353 


B19154 


-13 


48.26 


0.62 


94.5 


1.8 


5.7 


0.020 


0.017 


FW90 




J1922+2110 


B1920+21 


55.28 


2.93 


217.1 


4.8 


16.2 


0.045 


0.013 


this work 




J1926+1648 


B1924+16 


51.86 


0.06 


176.9 


5.2 


14.9 


0.034 


0.012 


this work 




J1930+1852 






54.10 


0.27 


308.0 


3.2 d 


10 d 


0.096 


0.031 


VB88; LAS01 


SNR:G54.1+ 0.3 


J1932+2020 


B19294 


-20 


55.57 


0.64 


211.2 


4.8 


14.9 


0.044 


0.014 


FW90 




J1932+2220 


B19304 


-22 


57.36 


1.55 


219.2 


10.4 


13.7 


0.021 


0.016 


FW90 




J1935+1616 


B19334 


-1(5 


52.44 


-2.09 


158.5 


5.2 




0.030 




FW90 




J1939+2134 


B1937+21 


57.51 


-0.29 


71.0 


4.6 


14.8 


0.015 


0.005 


FW90 




J1952+3252 


B1951+32 


68.77 


2.82 


45.0 


1.0 


4.0 


0.045 


0.011 


B84 


SNR:CTB80 


J2004+3137 


B20024 


-31 


69.01 


0.02 


234.8 


7.0 


12.0 


0.034 


0.020 


FW90 





a di and d M are lower and upper distance limits, respectively. A flat iFich et al.l (|l989l ) rotation curve with 
the standard IAU values of R a and o ( Kerr fe Lvnden-Bell 1986t) was used. See text for procedures to 
determine distances with other choices of R and o . 

b ( n e,u) and (n e ^) are upper and lower mean electron density limits, respectively. 

Reference kev: B84=|Blair et all dl984h: C75 ^Caswell et all dl975h: FW 90 =lFrail fc Weisberel (jl990h : 



LAS01= Euet all |200ll ; R97= lRosino et aD (jl997f ); VB88= IVelusamv fc Beckerl Jl988h. Other data are fro m 
the ATNF Pulsar Catalogue at http://www.atnf.csiro.au/research/pulsar/psrcat (jManchester et all 120051 ). 
Subsequent work bv llu et al. I (|2002h suggests d ~ 5 kpc under an assumption of uniform HI density in this 
direction, which the present authors find unlikely. (See Fig. 5.) 
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Table 3: Observing Parameters, Estimated Distances, and Transverse Speeds of High- 
Latitude Pulsars 



PSR J 


PSR B 


^obs 


-^~sys,off — line 


CT r,off — line 


DM 


d a 




Prop. Mo. 


Vtrans 


Vtrans 






(far) 


(K) 




(pc cm -3 ) 


(kpc) 


(mas/yr) 


Ref. c 


(AU/yr) 


(km/s) 


J0304+1932 


B0301+19 


1.3 


41 


0.015 


15.737 


0.62 


37 (5) 


LAS82 


23 


110 


J1239+2453 


B1237+25 


1.4 


41 


0.009 


9.242 


0.85 


115.6 (12) 


B03 


98 


470 


J1537+1155 


B1534+12 


1.7 


42 


0.13 


11.61436 


1.02 


25.086 (20) 


KWS03 


26 


121 


J1543+0929 


B1541+09 


2.9 


42 


0.03 


35.24 


3.49 


8.3 (10) 


B03 


29 


1 10 


J2305+3100 


B2303+30 


1.2 


41 


0.07 


49.544 


3.66 


20 (2) 


B03 


73 


350 



"Distance d is from the DM and the 



_ - .Cordes fc Lazid (2002) electron density model, except for J1239+2453 

[interferome t ric pa rallax; Brisken et al.l ( 2002h ] and J1537+1155 [differential galactic acceleration (see text); 
Stairs et all |2002h ]. 



fc The uncertainty in units of the last digit of the quoted value is given in parentheses. 
c Reference 



key: 



B03= 



Brisken et al 
1 — 



T 



LAS82= lLvne. Anderson fc Salter! (| 19821 ) . 

|http: / / www.atnf.csiro.au/research/pulsar/psrcat ( Manchester et al. 20051) 



(120031) : KWS03= lKonacki. Wolszczan. fc Stairs! (|2003l ): 
Ot her data are from th e ATNF Pulsar Catalogue at 
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Velocity (km/s) 

Fig. 1.— For PSR J1909+0254 = B1907+02. Top: HI emission spectrum in the direction 
of the pulsar. Middle: Pulsar absorption spectrum. The resolution is 1 km/s unless Hanning 
smoothing is indicated, in which case it is 2 km/s. Bottom: Radial velocity as a function of 



distance, calculated from a Fich et a" 



of R and O (jKerr fc Lynden-Bell 



(119891) flat rotation curve with the standard IAU values 



1986). See text for procedures to determine distances 



with other choices of R n and f 
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Fig. 2. — HI emission and pulsar absorption spectra, and rotation curve, toward PSR 
J1922+2110 = B1920+21. See Fig. 1 caption for additional details. 
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Fig. 3. — HI emission and pulsar absorption spectra, and rotation curve, toward PSR 
J1926+1648 = B1924+16. See Fig. 1 caption for additional details. 
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Fig. 4. — The mean electron density in the in ner galactic plane visible from Arecibo Ob- 
servatory. The IAU (IKerr fc Lynden-Belllll986l ) values of R a and O were used; see text for 
procedures to rescale with other values. Pulsars in Table 2 with di < 1 kpc or \b\ > 3° are not 
shown, (a) Measured electron density limits (error bars) and NE2001 model electron den- 
sities (circles), versus longitude, (b) The ratio of measured electron density to the NE2001 
model value, versus longitude. 
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Fig. 5. — The galactic plane, showing spiral arms (ICordes fc Lazioll2002l ) and Arecibo inner 
galactic plane pulsars with measured distances [upper and lower limits are delineated by the 
ends of lines pointing at the Sun, which is at (X,Y) = (0.0,8.5) kpc]. Galactic longitudes 
in the Arecibo range, 30° < / < 70°, are shown at the bottom and right edges of the plot. 
Note especially the long interarm path in the direction of PSR B1913+16, at / ~ 50°. 
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d (kpc) 




V r (km/sec) 



<5(=d/R Q ) 




V r (km/sec) 



Fig. 6. — Calculated radial velocity as a function of distance along a line of sight in the 
galactic plane at / = 50°. A flat rotation curve with the IAU (IKerr fc Lynden-Belll Il986l ) 
values of R Q and G Q yields the solid line, while a flat rotation curve with new galactic 
constants (see text) is shown by a dashed line, (a) Radial velocity versus distance d along 
the line of sight, (b) Radial velocity versus normalized distance 5 along the line of sight, 
where 5 = d/R . 
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Fig. 7. — HI emission and pulsar absorption spectra, and ro tation curve, towar d PSR 
J0304+1932 = B0301+19. The 1976-77 absorption spectrum of bickev etalJ Jl98lh is also 
plotted atop ours as a dotted line. Note that the absorption spectrum vertical scale is 
magnified for clarity. See Fig. 1 caption for additional details. 
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Fig. 8. — Changes in the absorption spectrum of PSR J0304+1932 = B0301+19 over twenty- 
two years, (a) The old and new spectra are shown after the new spectrum has been inter- 
polated to the same frequencies as the old spectra. This plot is similar to the central panel 
of Fig. [7J except that the new spectrum in the earlier figure was not resampled. (b) The 
difference spectrum, with a ±2<r error envelope superposed (dashed lines). 
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Velocity (km/s) 

Fig. 9. — HI emission and pulsar absorption spectra, and rotation curve, toward PSR 
J1239+2453 = B1237+25. See Fig. 1 caption for details. 
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Velocity (km/s) 



Fig. 10. — HI emission and pulsar absorption spectra, and rotation curve, toward PSR 
J1537+1155 = B1534+12. See Fig. 1 caption for details. 
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Fig. 11. — HI emission and pulsar absorption spectra, and rotation curve, toward PSR 
J1543+0929 = B1541+09. See Fig. 1 caption for details. 
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Fig. 12. — HI emission and pulsar absorption spectra, and rotation curve, toward PSR 
J2305+3100 = B2303+30. See Fig. 1 caption for details. 



